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A sample of 7253 n^e^e^ decay candidates with 1.0% background contamination has 

been collected by the NA48/2 experiment at the CERN SPS, which allowed a precise measure- 
ment of the form factor, the branching ratio, and the CP violating asymmetry of and K decay 
widths was investigated. The results of the — > n^e^e^ analysis, as well as the status of the 
K"^ ;r*/x+/i^ analysis based in the same data set, are reported. 
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Introduction 

The FCNC processes n^l^l^ (/ = e,iJ.) are induced at one-loop level in the Standard 

Model. Their decay rates are dominated by the long-distance contribution via one-photon ex- 
change, and have been described by the Chiral Perturbation Theory (ChPT). Several models of the 
vector form factor characterizing the dilepton invariant mass spectrum and the decay rate have been 
proposed [||, ^ The first observation of the —>■ Ti+e+e^ process was made at CERN more 
than 30 years ago followed by BNL E777 ^ and E865 ^ measurements. The most precise of 
these, E865, based on a sample of 10300 candidates, allowed a detailed analysis of the decay form 
factor and rate, and a test of the next-to-leading order ChPT calculation [|l|]. 

A recent precise measurement of the — > n^e^e^ decay based on the full data set collected 
in by the NA48/2 experiment at the CERN SPS and pubhshed in 2009 [0] is reported here. The 
status and prospects of the — > n^pL^pL^ analysis based on the same data set are also discussed. 

1. The NA48/2 experiment at CERN 

The NA48/2 experiment which took data in 2003-04, designed for charge asymmetry mea- 
surements, uses simultaneous and K beams produced by primary SPS protons impinging on a 
beryllium target. Charged particles with momentum (60 ±3) GeV/c are selected by an achromatic 
system of four dipole magnets with zero total deflection, which splits the two beams in the vertical 
plane and then recombines them on a common axis. The decay volume housed in a 1 14 m long 
cylindrical vacuum tank. Both beams follow the same path in the decay volume: their axes coincide 
within 1 mm, while the transverse size of the beams is about 1 cm. With 7 x 10^' protons incident 
on the target per SPS spill of 4.8 s duration, the positive (negative) beam flux at the entrance of the 
decay volume is 3.8 x 10'^ (2.6 x 10"^) particles per pulse, of which 5.7% (4.9%) are K+ {R-). The 
K+/K- flux ratio is 1.79. 

A description of the NA48 detector and 2003-04 data taking can be found in [^]. The decay 
volume is followed by a magnetic spectrometer housed in a tank filled with helium at atmospheric 
pressure, separated from the vacuum tank by a thin (0.31%Xo) Kevlar window. The spectrometer 
consists of four drift chambers (DCHs), two located upstream and two downstream of a dipole 
magnet which provides a horizontal momentum kick of A;? = 120 MeV/c for charged particles. 
The nominal spectrometer momentum resolution is Op/p = ( 1.02© 0.044 •;?)% {p in GeV/c). 

The spectrometer is followed by a plastic scintillator hodoscope (HOD) used to produce fast 
trigger signals and to provide precise time measurements of charged particles. The HOD is fol- 
lowed by a liquid krypton electromagnetic calorimeter (LKr) used for photon detection and particle 
identification. It is an almost homogeneous ionization chamber, 21Xq deep, segmented transver- 
sally into 13248 cells 2x2 cm^ each, with no longitudinal segmentation. Detectors located further 
downstream (hadron calorimeter, muon counter) are not used in the K"^ — > n^e^e^ analysis. 

2. n^e^e^ decay analysis 

The n^e^e^ rate is measured relative to — s- n^^n^ normalisation channel (where 71^ — > 

e^e^y is the Dalitz decay). The signal and normalisation final states contain identical sets of 
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charged particles. Thus particle identification efficiencies, potentially a significant source of sys- 
tematic uncertainties, cancel at first order. Three-track vertices (compatible with K"^ Ti^e+e^ 
and —>■ ti^ti^^ topology) are reconstructed by extrapolation of track segments from the spec- 
trometer into the decay volume, accounting for stray magnetic fields and multiple scattering. 

A large part of the selection is common to signal and normahsation modes, and requires a 
presence of a vertex satisfying the following criteria. 

- Vertex longitudinal position is inside fiducial decay volume: Zyertex > -Zfinai collimator- 

- The tracks should be in DCH, HOD and LKr geometric acceptance, and have momenta in the 
range 5 GeV/c < p <50 GeV/c. Track separations should exceed 2 cm in DCHl plane to suppress 
7 conversions, and 15 cm in LKr front plane to minimize effects of shower overlaps. 

- Total charge of the three tracks: Q = ±1. 

- Particle identification is performed using the ratio E / p of energy deposition in the LKr calorime- 
ter to momentum measured by the DCHs. The vertex is required to be composed of one n candidate 
{E/p < 0.85), and a pair of oppositely charged candidates (E / p > 0.95). 

If several vertices satisfy the above conditions, the one with the best vertex fit quaUty is con- 
sidered. The n^e^e^ candidates are selected by applying the following criteria. 

- Ti^e^e^ momentum within the beam nominal range: 54 GeV/c < \picee \ < 66 GeV/c. 

- n^e^e^ transverse momentum with respect to the beam trajectory (which is precisely measured 
using the the concurrently acquired — > 3?!* sample): pj < 0.5 x 10^^ (GeV/c)^. 

- Kinematic suppression of the main background channel — > n^n^ (and other backgrounds 
induced by tt^ and tt^^, 4e^ decays) by requiring z. = (Mee/MK)^ > 0.08, which corresponds to 
Mee > 140 MeV/c^, and leads to a loss of ~ 30% of the signal sample. 

- 7r*e+e^ invariant mass: 470 MeV/c^ < Mj^ee < 505 MeV/c^. The lower limit corresponds to a 
Ey < 23.1 MeV cutoff for the energy of a single directly undetectable soft IB photon. 

For the n^n^^ normalisation mode candidates, a presence of a LKr energy deposition 

cluster (photon candidate) satisfying the following principal criteria is required. 

- Reconstructed cluster energy E > 3 GeV, cluster time consistent with the vertex time, sufficient 
transverse separations from track impact points at the LKr plane. 

- e^e^ 7 invariant mass compatible with a decay: {Meey — Mj^o \ < 10 MeV/c^. 

- 7r*e+e^7total and transverse momenta: same requirements as used for —>■ n^e^e^. 

- 7r''=e+e" 7 invariant mass: 475 MeV/c^ < Mj^eey < 510 MeV/c^. 

The reconstructed n^e^e^ invariant mass spectrum is presented in Fig. [l]a. The number of 
K^e+e- candidates in the signal region is N^ee = 7253, of which 4613 (2640) are K+ {K^} 
candidates. The background sources are —>■ n^n^ and K"^ — > Ti^^e^v decays with tt^ —>■ e^e^j 
and jlt^ misidentification, and kaon decays with two or more e+e^ pairs in the final state from 
71° decays or external 7 conversions. Background contamination is measured to be (1.0 ±0.1)% 
using the sum of spectra of the unphysical data LFV %~^e^e^ and triple charge %^e^e^ candidates. 

The number of —>■ n'^Tip candidates in the signal region is N2n = 1.212 X 10'. The only 
significant background source is the —>■ 7r^/x*v decay, with contamination in the signal region 
estimated to be 0.15% by MC simulation. 

The decay is supposed to proceed through single virtual photon exchange, resulting in a spec- 



3 



Measurement of the K^l^l decays by NA48/2 



Evgueni Goudzovski 



























(a) 


































































r 



















£600 
> 

UJ 

500 
400 
300 
200 
100 

















(b) 




T 




































+ 

f 




































* 


















» 



" 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 

M(7t±e*e ), GeV/c' z 
Figure 1: (a) Reconstructed spectrum of n'^e^e^ invariant mass: data (dots) and MC simulation (filled 
area). Note the description of the radiative mass tail by the PHOTOS simulation, (b) z spectrum of the 
selected — > K^e^e^ candidates. Filled area: estimated background multiplied by a factor of 5. 



trum of the z = {Mee/MxY kinematic variable sensitive to the form factor W{z) [|I]]: 



dV 



dz 1271(471 



AV2(l,z,4)Jl-4^|l 



\W(z) 



(2.1) 



where = mg/Mn, rj^ = niji/MK, and A {a, b, c) = + b'^ + — lab — lac — Ibc. The following 
parameterizations of the form factor W{z) are considered in the present analysis. 

1. Linear: W{z) = G pM^f q{1 + 5z) with free normalisation and slope (/o,5). Decay rate and 
z spectrum are sensitive to |/o|, not to its sign. 

2. Next-to-leading order ChPT [[I]]: W{z) = GFM'j^{a++b+z)+W^^{z) with free parameters 
{a+,b+) and an explicitly calculated pion loop term W^'^(z) given in [|l|]. 



3. Combined framework of ChPT and \wge-Nc QCD 
W{z) = W(w,j8,z) with free parameters (w,j8). 



the form factor is parameterized as 



4. ChPT parameterization [^] involving meson form factors: W{z) = W(Ma,Mp,z). The reso- 
nance masses {Ma, Mp) are treated as free parameters in the present analysis. 



The Coulomb factor is taken into account following for instance 



. Radiative corrections to 
are evaluated with a PHOTOS [llO] simulation of the ^ n^f 



decay, and cross-checked with a generalized computation for a multi-body meson decay [Q]. They 
are crucial for the extrapolation of the branching ratio from the limited M^nee (equivalently, Ey) 
signal region to the full kinematic region: about 6% of the total — > ll^e^e~{y) decay rate fall 
outside the signal region Ey <l'h.\ MeV. 

The z spectrum of the data events (in the visible region z > 0.08) presented in Fig. [T]b. The 
values of dTj^ee/dz in the centre of each /-bin of z are computed as 



{dYnee/dz) 



Ni-Nf A2^{\-e2n) 1 n 



A;(l-£,) Az Xk 



BR(/^± n^n^) ■ BR(7r^). (2.2) 
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Figure 2: (a) dr„ee/dz (background subtracted, corrected for trigger efficiency) and fit results according 
to the four considered models, (b) Reconstructed spectrum of n^jx^jj.^ invariant mass: data (dots), 
7r*/i+/X^ MC simulation and ^37:^ background estimate (filled areas). 



Here Nj and N? are numbers of K^e^e^ candidates and background events in the /-th bin, 

N2n is the number of n^n^ events (background subtracted), A/ and are geometrical accep- 

tance and trigger inefficiency in the /-th bin for the signal sample (computed by MC simulation), 
A2it and £271 are those for —>■ 71^ events, Az is the bin width set to 0.02. The external inputs 
are the kaon lifetime %, and normalisation branching ratios BR(A'^ tt^tt^), BR(7r^). 

The values of dVjcee/dz and results of the fits to the four models are presented in Fig. ^a. 
The model-independent branching ratio BRmi in the kinematic region z > 0.08 is computed by 
integration of dTj^ee/dz, and differs from each of the model-dependent BRs computed in the same 
z range by less than 0.01 x 10^^. The differences between model-dependent BRs come from the 
region z < 0.08, as seen in Fig. |2|a. 

Systematic uncertainties due to particle identification inefficiencies, imperfect MC descrip- 
tion of the beamline, background subtraction, trigger inefficiency, radiative corrections, and fitting 
method are considered. The external uncertainties related to limited relative precision (2.7%) of 
BR(7r£) are also taken into account. 

3. Results, discussion and prospects 



The measured model-independent BRmi(z > 0.08), and the parameters of the considered mod- 
els are presented in Table 2. The correlation coefficients and 68% confidence contours for model 
parameters are shown in |^. Each of the considered models provides a reasonable fit to the data (as 
indicated in Fig. ^a). The data are insufficient to distinguish between the models. The measured 
form factor slope 5 is in agreement with earlier measurements [0, 0, 11], and disagrees to the me- 



son dominance models [12] which predict lower slope values. The measured |/o|, a+ and are 
in agreement with the previous measurement [|6|] ; a+ is in agreement with a theoretical prediction 



-0.6 



+0.3 
-0.6 



0131]. The measured w, j8 are in fair agreement with an earlier measurement | 
The branching ratio in the full kinematic range, which includes a model-dependence uncer- 
tainty, is BR = (3.11 ib0.04stat. ±0.05syst. ±0.08ext. ±0.07modei) X 10^^; it agrees to earlier mea- 
surements [Q, ||, |6|]. The DCPV charge asymmetry of decay rates is measured for the first time: 
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Table 1: Model-independent BRn,i(z > 0.08), and fit results for the considered models. 



BRmi X 10'^ = 


2.28 


± 


0.03stat. 




0.04syst. 




0.06ext. 


2.28 




0.08 


l/ol = 


0.531 


± 


0.012stat. 




0.008syst. 




0.007ext. 


= 0.531 




0.016 


8 = 


2.32 


± 


0.15stat. 




0.09syst. 






2.32 


± 


0.18 


a+ = 


-0.578 


± 


0.012stat. 




0.008.syst. 




0.007ext. 


= -0.578 


± 


0.016 


b+ = 


-0.779 


± 


0.053stat. 




0.036.syst. 




O.OHext. 


= -0.779 


± 


0.066 


w = 


0.057 


± 


0.005stat. 




0.004.syst. 




O.OOlext. 


= 0.057 




0.007 


i8 = 


3.45 


± 


0.24stat. 




0.17syst. 




0.05ext. 


3.45 




0.30 


Ma/GdW/c^) = 


0.974 


± 


0.030stat. 




0.019syst. 




0.002ext. 


0.974 


± 


0.035 


Mp/(GeV/c2) = 


0.716 


± 


0.011 Stat. 




0.007.syst. 


± 


0.002ext. 


= 0.716 


± 


0.014 


A(^:±J = (BR+- 


BR )/(BR+ 


+ BR ) = 


(-2.2±1.5stat. 


±0. 


6syst.) X 10 


^, corresponding to an 



upper limit of \^{K^ee)\ < 2.1 X 10"^ 90% CL. However the achieved precision is far from the 
SM expectation |A(7S:4J| ~ 10"^ % and even the SUSY upper limit of |A(^^^J| ~ 10"^ Q. 



An analysis of the — > K^jx^ jx'^ decay based on the same data sample is in progress. A 
sample of ~ 3100 decay candidates with 3% background from the — > 371* is selected. The 
n'^jjL^jX'' invariant mass spectrum is presented in Fig. Unlike the K^e^e^ case, the full 

kinematic z range is accessible, and effects of radiative corrections are suppressed. In addition to 
spectrum, rate and CPV measurements, the first measurement of the forward-backward asymmetry. 



which can be enhanced with respect to n^e^e in both SM and MSSM [ 15], is performed. 
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